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Level spacings and periodic orbits 
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Starting from a semiclassical quantization condition based on the trace formula, we derive a 
periodic-orbit formula for the distribution of spacings of eigenvalues with k intermediate levels. 
Numerical tests verify the validity of this representation for the nearest-neighbor level spacing (k = 
0). In a second part, we present an asymptotic evaluation for large spacings, where consistency 
with random matrix theory is achieved for large k. We also discuss the relation with the method of 
Bogomolny and Keating [Phys. Rev. Lett. 77 (1996) 1472] for two-point correlations. 
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The statistical distribution of quantum energy levels 
is conjectured to reflect the chaoticity or integrability 
of the underlying classical dynamics [[j]J^]. For classi- 
cally chaotic systems one expects spectral distributions 
like those computed in random matrix theory (RMT) 
whereas for classically integrable systems the quan- 
tum levels appear to follow the distribution for a Pois- 
son process . Attempts to explain this correspondence 
are based on Gutzwiller's semiclassical trace formula || 
which bridges the gap between classical and quantum 
mechanics. The prime result in this direction is Berry's 
analysis of the spectral rigidity based on the so- 
called diagonal approximation and the classical sum rule 
of Hannay and Ozorio de Almeida ||. Starting from 
a semiclassical quantization condition ||f7|], Bogomolny 
and Keating Q were able to extract information on the 
two-point correlations going beyond the results of jl| . It 
is the aim of the present work to follow a similar path 
in order to obtain semiclassical information on the level 
spacing distributions. 

The Gutzwiller trace formula can be expressed as a 
periodic-orbit sum for the integrated density of states 
N(E):=^Q(E-E n ), 



N T (E) ~ N(E) 
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in the semiclassical limit, i.e. for h — > 0. N(E) denotes 
the mean part following from Weyl's law, whereas the 
fluctuations, N^(E), are given by a sum over all primitive 
periodic orbits 7 and their r-fold repetitions, the prime 
on the second sum indicating the exclusion of r = 0. 
The (primitive) action and period of a periodic orbit are 
denoted by S 7 (E) and T 7 (E), respectively. The explicit 
form of the amplitudes A ir can be found in || for chaotic 
and in Q for integrable systems. We have given a ver- 
sion truncated at periods |rT 7 | = T, which will be needed 
in the following. Instead of directly using (|l|) to express 
spectral functions, we will use an approximate spectrum 
E n (T) obtained from the condition 



Before investigating the statistical distribution we have 
to unfold the spectrum such that its mean density d(E) = 
dN(E)/dE is rescaled to unity. To this end, inves- 
tigating spectral correlations in an interval I(E; K) :— 
[E — hw,E + huj], we introduce the unfolded energies 
x n {T) := E n {T)d, d := d(E), see, e.g., f(J for details. 
Spacings of two unfolded energies with k intermediate 
levels are given by s n (k;T) = x n+k+1 (T) - x n (T) and 
integrated level spacing distributions are defined by 



I(k,s;T) 



- E 

Ni ^ 



Q(s - s n (k;T)) , (3) 



E n ei(E-,h) 



where Ni denotes the number of eigenvalues contained in 
I(E; H). The often used level spacing densities P(k, s; T) 
are the derivative of I(k, s; T) with respect to s. Since the 
condition s > s„(fc;T) can be rewritten as NT(E n (T) + 
s/d) — NT(E n (T)) > k, we can substitute the argument 
of the step-function. Upon replacing the sum by an inte- 
gral over I(E; K) with weight d(E') which, in turn, we can 
asymptotically substitute by dr(E') := dN T (E')/dE' we 
obtain 



I(k,s;T) 



e (n t (e' + - Nt(E') -k-l 



d T (E'} 



(4) 



N T {E n {T)) = n + \ 



(2) 



The brackets denote the energy average (. . .) := 
iku Ie-Iiuj ■ ■ ■ dE', due to which we may asymptotically 
replace dT(E')/d by 1. Since the semiclassical limit 
h —* now implies d — > 00 , we have s/d <C 1 , and expand- 
ing the mean part of Nt about E yields the periodic-orbit 
formula 

J(M;T)~ (e(s-k-l + N$(E' + l} -jvH (£;'))} 

(5) 

for the integrated level spacing distribution. A similar 
expression can be obtained for the level spacing distribu- 
tion P(k, s; T) itself by taking the derivative of (||). Here 
we again neglect d T (E' + s/d)/d because of the energy 
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average, i.e. 

P(k,s;T)~(5 (s-k-l + N$(E' + J) -iV|(£')))- 

(6) 

For the case k = 0, i.e. for nearest-neighbor level spac- 
ings, we will now test these formulae numerically for two 
toy models. Our first example is a classically integrable 
system, namely a rectangular quantum billiard with as- 
pect ratio 2/(l + v / 5) and Neumann boundary conditions. 
For scaling systems the semiclassical limit H — > can be 
replaced by the high energy limit E — > oo, with a suitable 
change in the energy average. For the plots we have taken 
into account energies in the interval E e [0, 4000] (cor- 
responding to 214 eigenvalues). Figure |l| shows the in- 
tegrated level spacing distribution for this system where 
for the cut-off time T we have chosen T = Th '■= 2-Khd 
(dashed-dotted line) and T = 50 Th (dashed line), re- 
spectively. When increasing the cut-off time T we ob- 
serve convergence towards the quantum mechanical re- 
sult (solid line) which in turn in the limit E — > oo con- 
verges to a Poissonian level spacing distribution (dotted 
line), cf. @ 
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FIG. 1. Integrated level spacing distribution from pe- 
riodic orbits (H) for a rectangular billiard with cut-off 
time 50 Th (dashed), and Th (dashed-dotted). The 
quantum mechanical level spacing distribution (solid) 
and the graph for a Poisson process (dotted) are shown 
for comparison. 

For a second test we have chosen the imaginary parts 
of the non-trivial zeros of the Ricmann zeta function. 
These serve as a model for quantum chaos, since their 
spectral correlations are well described by the Gaussian 
unitary ensemble (GUE), see, e.g., and references 
therein. The density of the Ricmann zeros is related to 
a sum over primes, see, e.g., Jl3[ ], in a similar way as 
the density of states of a quantum system is related to 
a sum over periodic orbits of the corresponding classi- 
cal system (^). However, for chaotic systems as well as 
for the Ricmann zeta function it has been observed that 



there might be an optimal cut-off for the periodic orbit 
sum, i.e. including arbitrarily long orbits doesn't nec- 
essarily improve the result. In the case of the Riemann 
zeta function this is known as the Riemann-Siegel cut-off 
JTJ which corresponds to T H /2 [[TJJ. 

In figure ^| the integrated level spacing distribution of 
the first 649 Ricmann zeros with imaginary parts between 
and 1000 (solid line) is shown as well as the analogue 
of (||) (dashed line) with a truncation of the sum over 
primes at the Riemann-Siegel cut-off. We observe a good 
agreement of both curves at large values of s were both 
approach the GUE result (dotted line) . 
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FIG. 2. Integrated level spacing distribution for the 
zeros of the Riemann zeta function (solid) and the spac- 
ing distribution calculated from prime numbers (dashed) , 
cf. (||), T = Th/2. The spacing distribution of the GUE 
(dotted) and for a Poisson process (dashed-dotted) are 
shown for comparison. 

However, for values of s below 1 the periodic-orbit for- 
mula differs from both, the GUE and the exact curve. 
This is due to the fact that the periodic orbit sum (fif), 
truncated at the order of the Hciscnbcrg time, cannot 
reproduce features below the scale of mean level spac- 
ing. Unlike in the integrable case, we cannot simply 
increase the cut-off since now the sum over primes be- 
haves like an asymptotically divergent series, i.e. we ob- 
serve best agreement when truncating the sum at the 
Riemann-Siegel cut-off. Analogously we will adopt the 
same point of view for periodic-orbit sums of classically 
chaotic systems by treating the respective trace formu- 
lae like asymptotically divergent series, choosing T of the 
order of the Heisenberg time. 

We remark that in numerical calculations it is neces- 
sary to use a smoothed version of the periodic orbit sum 
(|l|) . We have chosen to add a so-called remainder term as 
introduced by Aurich and Steiner ||], which in the case 
of the Ricmann zeta function is related to Guinand's for- 
mula Q. 

Now that we have checked the periodic orbit formula 
numerically, we attempt to evaluate P(k, s; T) asymptot- 
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ically for generic quantum systems with chaotic classical 
limit. A well established conjecture on global eigenvalue 
statistics 17| states that (after suitable normalization) 



in the semiclassical limit the value distribution of 
is given by a Gaussian with zero mean. In analogy to 
H we also assume Gaussian behavior for the difference 
N^(E' + s/d)-N^(E'). Although this assumption leads 
to the correct result in the case of two-point correlations 
P,[i8| , we remark that it corresponds to neglecting cross- 
correlations between at different arguments, which, 
especially for small s, will become important. Thus, in 
the limit s — > oo, s/d— > 0, the energy average in (g) can 
be performed approximately, yielding 



P(k,s;T) 



1 
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: exp 
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2a 2 



(7) 



where the (s-dependent) variance a 2 :— ((N^iE'+s/d)— 

Nt(E')) 2 ) of the Gaussian is still to be determined. This 
last step can be done by substituting the trace formula 
(|l|) for Nj,, employing an expansion in s/d (cf. fig] ) and 
making use of the diagonal approximation p,H , 



\rT^\<T 
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where g denotes the generic multiplicity of orbits which 
share the same action, cf. p|Jl0| . Using the classical sum 
rule [El we obtain 



-7T 2 , 
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(10) 



7 denoting Euler's constant. We remark that the sub- 
leading term of this expansion is not unique at this 
point, since it is influenced by corrections to the sum 
rule, which unfortunately are unknown, as well as by the 
choice of the cut-off time. Therefore, when now setting 
T = CT H = C2ttM we obtain 
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(logs 



(11) 



Here a is kept as a free parameter, which, however, will 
not influence the large-s asymptotics. Upon substituting 
( pT| ) into (0) we obtain a semiclassical formula for the 
level spacing distributions as s — > oo, 



P SC (M 



1 n/(2g) 
log s + a 



exp 



^ 2 (s-fc-l) 2 
2.9 (log s + a) 



(12) 



Since ( [T2| ) is mainly supported around s w k + 1 we may 
for a moment substitute logs by log(fc + 1) for large k. 



The formula obtained that way is consistent with a con- 
jecture for P(fc, s) from RMT [jTJj, see also [^0|, which can 
be adopted to generic quantum systems with chaotic clas- 
sical limit ]2l| . However for small k and in particular for 
nearest-neighbor level spacings ([l2]) fails to reproduce the 
RMT results. This has to be understood in the following 
sense: The techniques used are essentially restricted to 
two-point correlations. Since the mean distance between 
Xn+k+i and x n is given by k + 1, the P{k,s) can only 
be dominated by two-point correlations on scales where 
s w k + 1. Thus an asymptotic evaluation as s — > oo also 
implies large k. 

We will now investigate the relation of our result ( p^ ) 
with the method of Bogomolny and Keating for the two- 
point correlation function i?2(s), which can be repre- 
sented as the sum 



i? 2 (s) :=J2P(k,s)-l. 



(13) 



fc=i 



Substituting (|T^) for the level spacing distributions and 
applying the Poisson summation formula yields 

-R2.se (s) w ~ exp(-27T 2 j/ 2 cr 2 ) exp(27rii/s) 



fez\{o} 



1 + erf 



s — 1 + 27rij/(T 2 



(14) 



where a is given by (|11|). The leading terms as s — > 00 
derive from v = ±1 (as can be easily seen from the first 
exponential), i.e. 



-R2.se (s) » 2e 



-2 9Q co s (2tts) 

S 2 9 



(15) 



Comparing with the large-s asymptotics of two-point cor- 
relation functions of RMT we observe that we have ob- 
tained the leading oscillatory contribution, which can be 
expected from the method of Bogomolny and Keating 
[pLl8[, but are missing the term which corresponds to 
the diagonal approximation of the spectral form factor, 
cf . . The reason for this is that between (||) and 

(g) we have neglected the term d%{E')d%{E' + s/d)/d 2 , 
which there was consistent, since it can be checked that 
doing so does not change the large-s asymptotics of (|l^). 
However, in the sum ([O]) the diagonal approximation of 
this term yields the missing term of R.2(s). 

In this sense our result for the level spacing distribu- 
tions (for large k) is in leading order consistent with RMT 
two-point correlations. Thus, we can now compare jl5| ) 
to the respective results from RMT, in order to fix a. 
For systems without time-reversal symmetry (g = 1) we 
have to compare with the GUE result, see, e.g. p3j, thus 
obtaining a = log(27r). Analogously for time reversal 
invariant systems (g — 2) we obtain a — log(v / 2 7r) by 
comparing with GOE two-point correlations. Thus, with 
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the substitution log s w log(A: + 1) for large k we have ob- 
tained Gaussians for the P(k, s) with approximate vari- 
ances 



a 2 (k) 
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(log(fc + 1) + a) . 



We briefly remark that the reasoning easily carries over 
to the third universality class, the Gaussian symplectic 
ensemble (GSE). To this end we have to consider trace 
formulae for particles with half integer spin ]24|,^5| . For 
time reversal invariant systems one has to take into ac- 
count Kramers' degeneracy as in . Again using the 
Gaussian ansatz (^) and methods of |10| for the calcula- 
tion of er 2 wc obtain 



log s + a 



exp 



7T 2 (S - k - If 

log s + a 



(16) 



where the same procedure as above yields a = log 8. 

Summarizing, we have derived periodic orbit formu- 
lae for the level spacing distributions P(k, s) and I(k, s). 
Numerical tests show that using purely classical input 
even the nearest-neighbor spacing distribution can be 
obtained for both, classically integrable and chaotic sys- 
tems, where in the latter case we cannot go below the 
scale of mean level spacing. In a second part we have 
presented an asymptotic evaluation of the formulae for 
large spacings which yields good approximations, if we 
restrict ourselves to the case of large k, but does not 
apply to the nearest-neighbor level spacing. 
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